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Abstract 

We study the effect of WIMP annihilation on the temperature of a neutron star. We shall 
argue that the released energy due to WIMP annihilation inside the neutron stars, might affect 
the temperature of stars older than 10 million years, flattening out the temperature at ~ 10^ K 
for a typical neutron star. 
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I. INTRODUCTION 



Since Zwicky proposed the problem of the "missing mass" in 1933, a lot of theoretical 
and experimental effort has been made in order to unveil the nature of dark matter. Today, 
WMAP has provided very accurate data regarding the matter density in the universe [1|. 
The energy density of the universe is composed of 4% atoms and roughly 22% dark matter. 
Data from recent observations indicate that dark matter cannot be attributed more than 
20% to dim objects like black holes, brown dwarfs and giant planets [^J- From a theoretical 
point of view, several candidates rise from different theories, such as neutralinos 



B3, 



Majorana neutrinos, and lately technibaryons provided by theories that are not ruled out 



by the electroweak precision measurements 



10|. From the experimental point 



of view, the focus is on the direct and indirect detection of dark matter particles. The direct 
detection might occur in underground experiments like CDMS that in principle can detect 
recoil energies from collisions between Weakly Interacting Massive Particles (WIMPs) and 
nuclei, or atmospheric experiments like XQC, where strongly interacting particles might 
collide with the detector. The indirect detection might occur via gamma-ray and neutrino 
telescopes, where the presence of WIMPs can be detected indirectly, by observing products 
of WIMP annihilations. In particular, provided that WIMPs can annihilate and because 
they can be trapped inside the earth or the sun, such annihilations would produce jets of 
particles and more specifically neutrinos coming straight from the center of the earth or the 
sun, that could be possibly detected by neutrino telescopes [Ul, ll2|, ll3| . On the other hand, 
gamma-ray telescopes can inprinciple detect gamma-rays produced by WIMP annihilation 
at the center of the galaxy ij] . Both direct and indirect detection experiments can impose 
strong constraints on the cross section of the WIMP with the nuclei. For instance, heavy 
Dirac neutrinos have been excluded as WIMPs for masses up to several TeV, because their 
elastic cross section with nuclei is sufficiently large and therefore they should have been 
detected by now in CDMS [l5|. 

In this paper we investigate the possibility of a different kind of indirect signature of 
WIMP annihilation. Instead of looking at the indirect signals from the annihilation of 
trapped WIMPs inside the earth or the sun, we examine the consequences of WIMP annihi- 
lation on the temperature of neutron stars. The neutron stars are massive compact objects 
with very low temperatures. Naively one might expect that since the mass of the trapped 
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WIMPs inside a neutron star represents a tiny fraction of the overall mass of the star, such 
an effect should be negligible. However, the annihilation of massive particles inside the star 
releases a huge amount of energy that is heating up the star. As we shall argue, once the 
accretion rate of dark matter particles equilibrates the rate of annihilation, the amount of 
released energy is independent of the star's temperature and therefore at late times the 
WIMP annihilation can keep the star at a constant temperature that depends on the mass 
and the radius of the star, the cross section of annihilation and the local dark matter density 
of the star. 

The paper is organized as follows: First we calculate the rate of dark matter accretion onto 
the neutron star including general relativity corrections. Then we calculate the annihilation 
rate for the WIMPs and we calculate the effect of the WIMP annihilation on the cooling 
curves of a typical neutron star made of regular nuclear matter. We present our conclusions 
in the last section. 



II. WIMP'S ACCRETION RATE ONTO THE NEUTRON STAR 



The accretion of dark matter particles inside the earth and the sun is not a new subject. 



Press and Spergel studied first in Uj the capture rate of WIMPs inside the earth and the 
sun. More elaborate calculations were also done by Gould 13] > taking into account 
several effects specifically for the case of the earth and the sun. An estimate of the accretion 
rate onto a neutron star was also provided by Goldman and Nussinov [l6|, who were the 
first to study effects of WIMPs on neutron stars. In this section we calculate the accretion 
rate of WIMPs onto a typical neutron star including also general relativity corrections that 



turn out to affect up to 70% the rate. Our derivation is along the lines of We assume 



that the WIMP population has a Maxwell-Boltzmann distribution of velocities 

3 . (-zv^ 



p{v)dv = no J Attv exp j dv, (1) 

where v = 270km/s, and uq is the number density of the WIMPs in the neighborhood of 
the neutron star. The flux of WIMPs (number per area per time) that crosses a spherical 
surface of radius R with velocity between v and v + dv and angle with respect to the normal 
between 6 and 6 + dO, is 

/ 3 /_3^2\ 



We can express the flux in a more convenient way with respect to the two invariants of the 
motion, i.e. the energy of the WIMP per mass E = {l/2)v'^ and the angular momentum per 



mass J = vRsinO. The total accretion rate (number of particles per time) is 11| 

dj' = AirR^dF = no( ^ exp ( Air^dEdJ^ (3) 

The actual capture rate of WIMPs by the star can be calculated in two steps. The first one 
is to determine what part of the phase space for E and J can give orbits for the WIMPs that 
intersect with the neutron star. In the second step we have to determine what fraction of 
the particles that intersect with the star, lose enough energy so they can be trapped inside 
the star. For the first part of the calculation, we have to find the trajectories that have a 
perihelion (closest distance to the center of the star) at most equal to the radius of the star. 
Press and Spergel calculated this using classical Newtonian mechanics. The perihelion is 



P'" \GM) ' \ V GMGM 

where G is the gravitational constant and M is the mass of the neutron star. This expression 
has two limiting cases. For J^E « (GAf)^ 

"■P^" " 2GM' 

and for J'^E » {GMf 



rperi = \/^. (6) 



The two regimes are separated by the hyperbola 

J2 E 



GM GM 



1. (7) 



Since Eq. (I3j) falls exponentially with respect to the energy, we approximate (as it is done 



in [HI]) the exponential as unity with E varying from zero to (l/3)'u^, which is the char- 
acteristic scale of the exponential. In addition, E is also restricted to values smaller than 
Eq, where Eq represents a constant that parametrizes the maximum kinetic energy per mass 
of the WIMP at asymptotically large distance from the star in order for the WIMP to be 



captured by the star. We shall determine Eq later on. Therefore, as it was argued in 



the accretion rate of capturable WIMPs is given by ([3]), if we integrate over E from zero 
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to the minimum between (l/3)t;^ and Eq and over from zero up to 2GMR (that comes 
from ([5]) if Tperi = R), where R is now the radius of the star. The rate can be written as 

:F=no (^^^ ^ A7t\2GMR) min Qti^, Eo^ . (8) 

This formula differs by a factor of 2 with respect to the corresponding one in llj, as it was 
first pointed out by Gould 12]. Although the above formula is a good approximate relation 
for the capture rate of WIMPs for the sun and the earth, for the case of a neutron star, 
general relativity corrections increase the rate significantly. To show this, we are going to 
use the timelike geodesic equations that describe the motion of a particle in a Schwarzschild 
metric. The trajectory for nonrelativistic particles (as the WIMPs) is given by 17] 



du\^ n ^ 2 2E 

-j =2m«'-« +^m« + -^, (9) 

where m = GM (in natural units) and u = 1/r. We want to find for what values of the phase 
space of E and J^, the perihelion becomes smaller or equal to the radius of the neutron star. 
At the perihelion du/d(j) = and u = 1/R (for rperi = R). If we express E in units of GM/R 
and in units of GMR, Eq. gives 

1 / 2GM\ , , , 

The above equation gives the relation between E and J in order the perihelion to be R. For 
E = 0, J"^ = 2/{l — 2GM/R). This means that the allowed phase space for has increased 
compared to the Newtonian case from 2 (in units again of GMR) to 2/(1 — 2GM/R). 
For a typical neutron star of mass 1.4 the solar mass Mq and a radius of 10 km, = 
2/(1 — 2GM/R) ^ 3.4. This a 70% increase in the phase space of and the capture rate 
compared to the classical case. Therefore Eq. should be modified as 

Now we estimate Eq. We shall show that in the case of a neutron star, £"0 >> (1/3)?;^ 
and therefore Eq. (ITT]) should be always taken with {l/3)v'^ as the minimum. A WIMP that 
intersects with the neutron star might or might not interact with the nuclear matter inside 
the star. If it does scatter at some point, the recoil energy and the energy loss of the particle 
is < T < 4m„m^/ (m„ + m-^)^, where T is the recoil energy and m„ and are the masses 



of the nucleus and the mass of the WIMP respectively. If we assume that the scattering is 
isotropic, then the recoil energy should be uniformly distributed. The condition that holds 
in order to capture a WIMP is that the energy loss in the scattering should be at least equal 
to the initial kinetic energy of the WIMP at asymptotic large distance from the star. If 
this condition is fulfilled, the WIMP stays in a bound state with the star. Therefore, for an 
average collision that takes place in the star, this condition can be written as 



where we took into account the gravitational redshift effect. Again we have chosen to set 
c = 1. If we plug the typical values we used before for the mass and the radius of a neutron 
star and a WIMP mass of the order of TeV, we find that Eq is three orders of magnitude 
larger than (l/3)i7^ and therefore for all the cases of interest, we can use ffTTl) with (l/3)-u^ 
as the minimum. 

Eq. fITT]) gives the rate of capturable WIMPs, that is the number of WIMPs per second 
that intersect with the neutron star. However, as mentioned before, in order for the WIMP 
to be trapped in the star, one or more collisions have to take place. We know from classical 
mechanics that if the WIMP does not scatter while travelling through the star, it cannot be 
captured by the star. We calculate now the fraction of the capturable WIMPs (given by (fTT!) ). 
that can yield scatters into bound orbits. Such a fraction would depend strongly on the 
elastic scattering cross section of the WIMP-nucleus system. For the typical neutron star of 
mass 1.4Mq and R = 10 km, the average density of neutrons is p = 3M/ (47ri?^m„) ~ 4x 10^* 
neutrons/ cm^. If we take a typical value for the elastic cross section between WIMP-neutron 
of the order of 10~^^ cm^, the mean free path is about 1 km. Since for an average WIMP, 
even one scattering is enough to result in a bound orbit around the star, a mean free path 
of 1 km means that the fraction of the capturable WIMPs that will be trapped is very close 
to 1. To entertain this, if for simplicity we assume that the WIMP has a straight trajectory 
while inside the neutron star, a segment of 1 km corresponds to an impact parameter of 9.9 
km, which means that only if the WIMP intersects between 9.9-10 km from the center of 
the star will travel a distance less than 1 km inside the star. Obviously if the cross section 
is larger than 10~^^ cm^, the fraction saturates even faster to 1. We shall give now a more 
quantitative answer about the fract.on followmg the der.vat.on of The fraction / of the 




(12) 



particles that undergo one or more scatterings while inside the star is defined as 



/ = ( 1 - exp 



m.. 



dl), (13) 



where the last approximation holds if the elastic cross section between WIMP-nucleus is 
smaller than o"crit = mnR'^/M ~ 6 x 10~^^ cm^. Eq. ( fT3l) now reads 

f^^ll-^-). (14) 



In order to find /, we examine trajectories with = as in [ll|], since E « GM/R. We 
average over (that ranges from to 3.4 for our typical neutron star). For an accurate 
calculation of /, we need to know the exact density profile of the star, in order to know 
explicitly the mass of the star M(r) as a function of the radius. Here we give an estimate 
of /, by assuming for simplicity that the density of the star is constant through the whole 
volume. This means that M{r)/M = (r/R)^ (where M is the total mass of the star). If we 
take the derivative of Eq. (Q with respect to 0, we get the following equation of motion for 
the WIMP inside the neutron star 

J2 M +^ i? M ~ J^u^^ Ru ' ^ 

where u = Ru and again is measured in units of GMR. The initial condition of this 
differential equation are u(0 = 0) = 1, which means that we have chosen </> = at the point 
where the WIMP crosses the surface of the star. For the velocity 



dui(l) = 0) _ I2GM ^^g^ 



R J2 

The length of the path of the particle travelling inside the star is 



dl 
R 




(17) 



We can find the length of the path of a particle inside the star if we integrate (j) from up 
to the angle that u{4>) = 1 again, which is the point where the particle exits from the star. 
Using Eqs. ( fT4l) . ( ITTI) . and after having solved Eq. ( |T5i) numerically, and averaging over 
from to 3.4, we found that the average path inside the star is 1.87i? and / = 0.45(T^/(Tcrit- 
We should emphasize that the above estimate of / holds for < o"crit- It is understood 
that if > Ucrit, / increases, saturating to 1 as soon as (Xcrit becomes larger than roughly 



10~^^ cm^. In addition, we should mention that our estimate is a lower bound for /. This 
is because in our derivation we assumed a constant density. Neutron stars are expected to 
be denser as we approach the center. In this case, / will be larger than our estimate. To 
illustrate this, we examined the extreme case where all the mass of the star is concentrated 
at the core. In this scenario, although somewhat unrealistic, / increases drastically because 
for a large range of J^, the orbits follow spirals around the center. A large fraction of the 
particles will be trapped in the star not because of energy loss due to a collision, but due to 
the fact that they are trapped gravitationally. Our expectation is that the real case should 
be somewhere in the middle and therefore we consider our previous estimate for / as a lower 
bound. 

Using Eq. (fTTj) and the values for our typical neutron star M = 1.4Mq, R = 10 km, we 
get the rate of accretion of dark matter inside the star in particles per second 

^ 3.042 x 10^5 ^ ^ 
m^(GeV) 

where A is a constant that parametrizes the local dark matter density in the vicinity of the 
neutron star in units of 0.3 GeV/cm^ (which is the standard dark matter density around the 
earth). For cross sections > 10"^^ cm^, / = 1, otherwise / is given by / = OAbay-facrit- 



III. ANNIHILATION RATE 



Once the WIMP undergoes one scattering inside the star, it loses on the average enough 
energy to be captured by the star. Even if the kinetic energy is sufficient enough to make 
it exit from the star, it will be forced to return and probably scatter again loosing even 
more energy. The WIMP can repeat this process several times until its kinetic energy 
reduces down to the thermal velocity inside the star. It is easy to show that for most 
cases of interest, the WIMP thermalizes very fast compared to the other time scales of the 
problem. We can make a very rough estimate of how long it takes for a WIMP moving 
with the average velocity of 270 km/sec to obtain a velocity comparable to the thermal 
velocity. Let's assume that the WIMP has undergone one scattering and therefore has lost 
on average energy of (2m„/m^)f^g(,/2, where ^esc is the escape velocity from the star. A 
simple approximate calculation shows that it will take a few seconds before the WIMP 
intersects again with the star and loose again a fraction 2m„/m^ of its energy. This depends 
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on the cross section of elastic collision between WIMP-nucleus, but as we showed in the 
previous section, for cx^ > 10~^^ cm^, this will happen on average. In that case, it will 
take just a few hours before the kinetic energy of the WIMP reduces down to the thermal 
velocity. So our conclusion is that for not extremely small cross sections, captured WIMPs 
will thermalize pretty fast and they will have a Maxwell-Boltzmann distribution in velocity 
and in distance from the center of the star. 

The population of dark matter WIMPs inside the star is governed in principle by three 
processes. The first one is the accretion of WIMPs onto the star. The second is the evap- 
oration and the third one is the annihilation. Once the WIMPs thermalize inside the star, 
they follow a Maxwell-Boltzmann distribution in velocity. Particles that are in the tail of 
the distribution (with large velocities) can escape from the star, if the velocity is larger 
than the escape velocity of the star. However, in the cases we are interested in, this 
process is exponentially suppressed. The rate of evaporating particles is proportional to 



exp{—GMm^/RT) [18|. Since the radius of a neutron star is very small, and we are in- 
terested in WIMPs with mass of the order of TeV, the rate becomes negligible. For a 
temperature of 100 keV (which corresponds to a typical temperature of a neutron star that 
is a few thousand years old), and for = 1 TeV, the suppression is exp(— 10^). Therefore 
we can safely ignore the evaporation process. 

If the WIMP is a Majorana particle, for example a Majorana neutrino, it is possible to 
co-annihilate with another one. The annihilation depends on the cross section as well as the 
density of the WIMPs inside the star. The annihilation cross section should not be confused 
with the elastic cross section between WIMP-nucleus that was mentioned before. One big 
difference between the two is that the annihilation cross section for Majorana particles is 
usually velocity dependent. If the WIMP is a Majorana neutrino, it has an elastic cross 
section with nuclei 

a, = ^^^/., (19) 

TT 

where /i is the reduced mass of the system WIMP-nuclei and Is is a form factor that depends 
on the nuclei [l^. The annihilation cross section of two Majorana neutrinos depends on what 
channels are open for annihilation. If the mass of the Majorana neutrinos is larger than 100 
GeV, the dominant channel is the annihilation to a pair of — W~ mediated by a Z 



boson [20|. In the case where >> 100 GeV, the annihilation cross section is given by 

671 

where /3 is the velocity of the WIMP (at the center of mass frame). Once the WIMP will 
get thermalized in the star, = 3T/(2m^) (T being the temperature in the star). 
The annihilation rate of WIMPs in the star is given by 

= (a^v) I n^dV, (21) 



where (cx^f ) is the thermally averaged annihilation cross section times the velocity and n is 
the density of the WIMPs inside the neutron star. For convenience we shall assume that the 
density is constant inside the star. In this case, the population of WIMPs inside the star is 
governed by 

^=J^-CaN'. (22) 
at 

The constant Ca = (c"x^)/^' "^here V is the volume of the star. The accretion rate JF was 
derived in the previous section. The solution of Eq. (12^ is 



A^W = A/£tanh(^). (23) 



The time scale r = 1/ The released energy due to the annihilation is 

E = CAN'^rn^ = Timi\^{tlT)m^. (24) 

The amount of the released energy depends on the time scale r. If r is large compared to the 
age of the known neutron stars, the hyperbolic tangent is suppressed and the effect of the 
dark matter on the temperature of the star is negligible. For particles that the annihilation 
cross section is velocity independent, r is given by 

2.1 X 10^ years 



(25) 



where (J39 is defined through (^oav) = usglO"^^ cm^. In the case of a Majorana parti- 
cle, Eq. fl2^ holds only approximately because the cross section is velocity dependent and 
therefore in thermal equilibrium temperature dependent. Since the temperature of the star 
changes as a function of time, this means that there is additional time dependence on the 
annihilation rate. Generally, for a Majorana neutrino with mass larger than 100 GeV, the 

in 



cross section is given by Eq. (!20|) times a factor sin^ 6 which denotes the suppression of the 
cross section due to a mixing between the left handed neutrino (that interacts weakly) and 
a sterile right handed neutrino p^, 120(. The time scale r is 

_ 1 _ 5.98 X lO^years , . 

\ ml KlQi) 



where the temperature T is measured in Kelvin degrees and (T39 is defined by the relation 
o"A = cTsglO"^^/?^ cm^. The mass is measured in GeV. Eq. fl2B|) can be written more 
conveniently in terms of the mixing angle sin 9 as 

2.52 X lOVears ^^7) 

If the Majorana neutrino is exclusively left handed, sin 6' = 1 and for a temperature of 10^ 
Kelvin, the time scale is about 10^ years (depending on how much larger is the local dark 
matter density in the vicinity of the star compared to the one of the earth). As it can be 
seen from ( |2^ . the annihilation saturates to JF for times larger than roughly 3r. As we 
already mentioned, since the temperature of the star changes in time, Eqs. (126|) and (!27|l 
are approximate. We shall return to the question of how fast the annihilation rate reaches 
the saturated value in the next section. 



IV. COOLING AND HEATING THE NEUTRON STAR 

In this section we investigate the influence of the WIMP annihilation on the temperature 
of the neutron star. Naively, one would expect that such an effect should be negligible 
due to the fact that the accretion of dark matter represents a tiny fraction of the whole 
mass of the star. However, there are two elements that make this investigation interesting. 
The first one is that although the trapped dark matter represents a small fraction of the 
mass of the star, the annihilation of two WIMPs releases a huge amount of energy. After 
the annihilation, this energy is carried mostly by leptons, quarks and photons. Since they 
cannot escape from the star, they will heat it up. A small portion of the energy will be 
carried by neutrinos that will escape and they will not contribute to the heating of the star. 
However, to first approximation, the energy carried by the neutrinos is negligible compared 
to the one carried by quarks, leptons and photons. So, we are going to assume that the 
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whole annihilation energy will be carried not by neutrinos and therefore this energy will 
heat up the star. The second reason we investigate this effect is that the energy released 
by the annihilation and consequently the emissivity of this process does not scale with the 
temperature. As long as equilibrium between accretion and annihilation has been reached, 
the released energy per time remains unchanged. All the dominant cooling processes of 
a neutron star scale with positive powers of T. This means that inevitably, even if the 
emissivity due to WIMP annihilation is small, as the temperature of the star decreases, the 
WIMP annihilation emissivity will dominate at some point. 

Let's assume for the moment and we shall examine later under what conditions this is 
possible, that the time 3r has been reached and the released energy from the annihilation 
of the WIMPs is E = Tmy.. The emissivity, i.e. released energy per volume per time is 

There are several processes that contribute to the cooling of a neutron star depending on 
the form of matter, the temperature, and the density of the star. If the star is sufficiently 
dense, unpaired quark matter or other exotic phases might occur deeply inside the star. If 
the neutron star has unpaired quark matter, for the first million years cools very fast due to 
neutrino emission via the direct Urea process. In this case the emissivity scales as ej, ~ T^. 
Direct Urea processes are allowed in sufficiently dense nuclear matter, nuclear matter with 
pion condensation, kaon condensation, or nonzero hyperon density, and in all phases of quark 



matter except CFL (2]| and references therein. For neutron stars that are not in sufficiently 
high density, direct Urea processes n— ^p + e + i>, p + e— s>n + i/ are kinematically forbidden. 
In this bystander neutron is needed in order to assist kinematically the reaction. 

This is the so-called modified Urea process. During the epoch dominated by the modified 
Urea, the star loses energy through neutrino emission, by converting protons and electrons 



22| 



to neutrons and vice versa. The emissivity of this process scales as ~ T*^. It is 

e.Hl.2xl0^ergcm-3s-)y j , (29) 

where n is the baryon density of the star and tiq = 0.17 fm~'^ is the baryon density in nuclear 
matter. In our calculation for the neutron star of M = IAMq and R = 10 km, we are going 
to use the average density n = 3.9 x 10'^'^ particles per cm^ and therefore n/uQ = 2.3. 

After the first million years and roughly as soon as the temperature of the star drops 
below 10^ K, the dominant mechanism of cooling is not anymore through neutrino emission, 

19 



but through photon emission from the surface of the star. The rate of heat loss from the 
surface of the star is 

L, = 47ri?VTlface, (30) 

where a is the Stefan-Boltzmann constant and Tgurface is the temperature of the surface of 
the star. The surface of a neutron star is usually colder than the interior of the star. This 
change in the temperature occurs inside the crust of the neutron star, taking place within 
100 meters below the surface. The surface temperature of the star is well approximated 



by [23|,y,l25| 



r.„,„ = (o.87><io»K)(^^^j (lOTf) ■ (31) 

where T is the interior temperature of the star and Qs = GM/ is the surface gravity. The 
rate of heat loss can now be expressed in terms of the interior temperature as 

L,^4.«V(0.87xlO»Kr(^)(^)'^ (32) 

If we divide over the volume of the star, we can get an "effective" emissivity of photons 
measured in energy over volume and time 



2.2 



'^ = m^'-'''''° [Wk) «Scm-=-'. (33) 

where we used gs = 1.85 x 10^^ cm/s^. 

In order to be able to derive the temperature as a function of time, we need to know the 
heat capacity of the star. For a gas of noninteracting fermions, the specific heat is given 



by |22| 



= E + fe)^ (34) 

i 

where the sum runs over the different species. In the case we investigate, namely the one of 
noninteracting nuclear matter, i runs over n, p, e and the Fermi momenta for neutral matter 
in weak equilibrium are 

Pi = (340 MeV) l—j (35) 

=PF = (60MeV)(— j . (36) 
The cooling of the star is dictated by the differential equation 
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Internal Temperature 
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FIG. 1: Left Panel: The internal temperature of a neutron star (in Kelvin) with M = IAMq and 
R = 10 km as a function of time (in years). The solid line that crosses the time axis corresponds to 
the case where the effect of dark matter annihilation is neglected. The thin solid line corresponds 
to a local dark matter density for the star of 0.3 GeV/cm^. The dashed and the thick solid lines 
correspond to local densities of 3 and 30 GeV/cm'^ respectively. Right Panel: As in the left panel 
for the surface temperature of the neutron star. 

dT ^ —Li, — + Ldm _ V{—ey — e-y + Cdm) _ — + edm 

Vc^ ~ Vc^ ~ ' ^ ' 

where the volume of the star V drops out at the end. We have neglected the contribution 
of the WIMPs to the specific heat, since they consist a tiny fraction of the mass of the star. 
We solved Eq. ( !37l) numerically, by imposing an initial temperature for the star of 10^° K at 
very early time. However, we should emphasize that the temperature is very insensitive to 
the initial condition. It affects only the temperature during the first years of the star's life, 
but it has no effect later on. In Fig. 1, we have plotted in a logarithmic scale the internal 
and the surface temperature of the star as a function of time, starting from time t = 1000 
years, up to 100 million years. We have plotted the temperature for 3 different cases that 
correspond to different local dark matter densities for the vicinity of the star. We chose the 
local dark matter densities to be the one of the earth (0.3 GeV/cm^), 10 times larger and 
100 times larger the one of the earth. For comparison we also plotted the cooling curve of 
the same star without including the effect of dark matter annihilation. As we can see from 
the figures, the dark matter annihilation does not affect the temperature of the star up to 
t = 10 million years. Between 10^ and 10^ years, the star cools due to neutrino emission 
via the modified Urea process, while after roughly 1 milion years through photon emission 
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from the surface of the star. However, by inspection of Eqs. (128|) . (p9l) . and (133!) . we see that 
the dark matter annihilation emissivity scales with the lowest power of T. More precisely, 
Cdm does not depend on T, once equilibrium between the rate of accretion of dark matter 
and the rate of annihilation has been established. This means that inevitably, when the 
temperature of the star drops sufficiently, the power of the dark matter annihilation that 
heats up the star will equate the power of photon emission and as a result the temperature 
will remain flat as a function of time. This happens roughly at t = 10 million years and at 
surface temperatures between 3000 to 10000 K (depending on the local dark matter density, 
the mass and the radius of the star). To entertain the possibility of having a neutron star 
with a local dark matter density 10 or 100 times larger than 0.3 GeV/cm^, we can use a 
indicative proflle density for the dark matter halo. We consider the Navarro-Frenk- White 



profile, where the dark matter density is given by 
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Pir) = ^- as^. (38) 

This profile has a spike in the center of the galaxy for positive 7. We shall use a = 1, /3 = 3, 
7 = 1, i? = 20 kpc and po = 0.235 GeV/cm^. Given this density profile, a neutron star with 
density 10 and 100 times larger than the local dark matter density of the earth, should be 
1.37 and 0.15 kpc from the center of the galaxy respectively. The position of the earth is 
roughly 8 kpc from the center of the galaxy. This means that a neutron star that exhibits 
the fiatness in temperature of ~ 10^ K for time t > 10 million years due to the dark matter 
annihilation, should lie at least ~ 6.5 kpc away from the earth. This limit can be improved, 
if the star has "more convenient" mass and radius from what we have considered for a typical 
neutron star. The emissivity edm is proportional to JF, which is proportional to the factor 
MR/ (1 — 2GM/R). The emissivity edm is also proportional to the local dark matter density. 
When we quote results with local dark matter density 10 times larger the one of the earth, 
this does not imply per se that the star has to be at a region of the galaxy with pdm = 3 
GeV/cm^, but the factor MR/{1 — 2GM/R) times the local dark matter density should be 
ten times the same factor for our typical star with M = 1.4Mq and R = 10 km times the 
density of 0.3 GeV/cm'^. For example, a neutron star of M = 2Mq and R = Q km, gives a 
factor of ~ 26.5 compared to our typical star. Therefore our results for density 100 times the 
density of the earth are applied also for a star with mass and radius given in the previous 
sentence and a local dark matter density of only 100/26.5 = 3.77 times the density of dark 

in 



matter around the earth. 

Now we return to answer a question we posed earher. In the results we have presented, 
we have assumed that equihbrium between the accretion and the annihilation rate has taken 
place before t = 10 million years, where the effect of the WIMP annihilation becomes im- 
portant. If r is much larger than 10 million years, the annihilation rate has not reached the 
asymptotic value ~ JF and the effect on the temperature of the star will be negligible. We 
saw in Fig. 1 that even if the asymptotic annihilation rate is reached very early, the WIMP 
annihilation does not affect the temperature for t < 10 million years (or equivalently for 
temperatures higher than 10^ K). This means that up to 10'' years, the temperature of the 
star is controlled by modified Urea and photon emission. We already mentioned that if the 
WIMP annihilation cross section is velocity dependent (like in the case of Majorana neutri- 
nos), Eqs. dMj), fl26|l . and fl27|) hold only approximately since r is temperature dependent 
(and implicitly time dependent). In particular, Eqs. (l26l) . and (1271) can give a lower bound 
estimate of how fast the asymptotic annihilation rate is reached. For most candidates of 



our interest, Majorana neutrinos, or Majorana technibaryons [10[, we can safely take / = 1, 
since the elastic cross section with nuclei is not smaller than 10~^^ cm^. By using Eq. (1271) . 
with / = 1, y4 = 10, and an "average" temperature of 10^ K for the star, we get that for cross 
section suppression sin 6' > 0.1, r ~ 10^ years. In reality, the situation is better, because 
this estimate is done assuming constant T = 10^ K. However, the temperature of the star is 
much higher at the beginning (~ 10^° K) and due to modified Urea drops down to T = 10^ K 
roughly at t = 10^ years. The higher temperature at the beginning, brings the annihilation 
rate faster to the asymptotic value. In addition, for 10^ < t < 10^ years, where the star 
cools mainly due to photon emission, the annihilation rate can also improve towards the 
asymptotic value, although since the temperature falls really fast, this happens with a slow 
rate. For the two cases we mentioned, Majorana neutrinos and Majorana technibaryons, 
the annihilation rate reaches the asymptotic value before 10 million years. For example for 
a majorana neutrino, Eq. ( 1271) . for sin 6* = 1, A = 100, and / = 1, gives 2.52 x 10^ years (for 
T = 10^ K). For a Majorana neutrino (or Majorana technibaryon) with suppressed coupling 
to the Z boson in order to account for the right dark matter density, for a mass of 1 TeV, 



sin^ = 0.26 



10|. In this case, Eq. (1271) gives a characteristic time scale r = 3.7 x 10^ years. 



with A, /, and T as before. Therefore, for most cases of interest and unless the cross section 
is very small, the annihilation rate reaches the accretion rate before 10 million years, which 



is the time where the dark matter annihilation affects the temperature of the neutron star. 



V. CONCLUSIONS 

We investigated the effect of WIMP annihilation on the temperature of a neutron star. 
We found that for a typical neutron star with a local dark matter density at least 3 GeV/cm^, 
and if the WIMP has an elastic cross section not smaller than 10~^^ cm^, the WIMP anni- 
hilation flattens out the temperature of the star around ~ 10^ K at t = 10 million years. 
Two neutron stars that have different local dark matter densities would have different final 
temperatures scaling as T ~ Pdm^- Given the uncertainty in our knowledge of the age of 
a neutron star and the fact that the peak of a blackbody spectrum of 10^ K lies in the 
infrared, it is a challenge to observe such an effect, which would possibly be a signature 
of WIMP annihilation. Alternatively, instead of trying to spot neutron stars with such a 
low temperature, it might be more efficient to study pulsars, detected already by their non- 
thermal emission and to constraint their thermal emission putting an upper bound on their 
temperature. 

We disregarded in our analysis the interesting possibility of having an exotic phase of 
quark matter in the typical neutron star. Exotic quark matter phases like color supercon- 
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30|. This could 



ductivity can have an effect on the cooling of a neutron star 
change the dominant cooling process for the first million years, but it is highly unlikely 
that our conclusions would change regarding the effect of the WIMP annihilation. We also 
neglected reheating mechanisms that might possibly be present for old neutron stars. These 
mechanisms can be viscous dissipation of rotational energy within the star {si], energy re- 



lease due to weak deviations from beta equilibrium 32], accretion from interstellar gas, or 



other [33| and references therein. All of these mechanisms are model dependent and it is not 
clear what is their effect on the temperature of an old neutron star compared to the WIMP 
annihilation mechanism. As a rule, in a nonsuperfiuid old star without a magnetic field and 
an accreted envelope, the alternative mechanisms mentioned above cannot probably prevent 
the star from getting a very low temperature. In this case, WIMP annihilation might be 
the dominant reheating mechanism of the star. 

Much theoretical work remains to be done. An interesting question is to study if the 
effect of WIMP annihilation inside a neutron star can change the suggested mechanisms 
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for superbursts and gravitational wave bursts proposed in [3J, |35|, |36| . These mechanisms 
assume a phase transition to quark matter at the center of the neutron star via metastable 
phases, that can lead to a burst. The released energy from the WIMP annihilation inside 
the star might work as a catalyst, accelerating the phase transition, making the metastable 
phases impossible to persist. 
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